Triplet superconductivity in the skutterudite PrOs4Sbi2 
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Abstract 

There is mounting evidence for triplet superconductivity in the recently discovered skutterudite 
compound PrOs4Sbi2- In this work, we propose nodal order parameters for the A- and B-phases 
of this superconductor which are consistent with angle dependent magnetothermal conductivity 
measurements and with low-temperature thermal conductivity data in the range T > 150 mK. The 
quasiparticle density of states and the thermal conductivity k zz are derived within the quasiclassical 
approximation . 

PACS numbers: 



1. Introduction 

The skutterudite compound PrOs4Sbi2 is a heavy-fermion superconductor with a transi- 
tion temperature of 1.8K. Jl|,|2, |3| Angle dependent magnetothermal conductivity measure- 
ments on this material have revealed an interesting multi-phase structure, characterized by 
energy gap functions A(k) with point nodes. |^,|5j Previously, several unconventional order 
parameters, including s+g-wave symmetry, have been proposed to account for this nodal 
structure, p Most of these models that have so far been considered for PrOs4Sbi2 are not 
able to describe self-consistently the observed angle dependence. J7|] 

Recently, there has been mounting experimental evidence for triplet superconductivity 
in this compound. First, from /i-SR measurements Aoki et al. discovered a remnant mag- 
netization in the B-phase of this compound, indicating triplet pairing. 8( Second, thermal 
conductivity measurements along the [0 1] direction at low temperatures (T > 150 mK) 
indicate T-linear and H-linear behavior, consistent with a triplet order parameter. J9| Third, 
the observed angle dependence of k zz with H rotated within the x-z plane indicates triplet 
superconductivity. Finally, recently reported NMR data by Tou et al. 10] for the Knight 
shift in PrOs4Sbi2 also suggest triplet pairing. 

The objective of this paper is twofold. First, we propose a new set of nodal order pa- 



rameters for the A and B phases o: 
conductivity data observed in Refs 



PrOs4Sbi2, which are able to describe the thermal 



The proposed model can also describe the 

n 



i and li 

isotropic superffuid density in the B-phase reported by Chia et al.[U|, albeit with the provi- 
sion that it has to be assumed that the nodal points in this experiment are aligned parallel 



13] 



to the external magnetic field when the sample is field-cooled. |12L Il3| In order to interpret 
the observed ^-dependence of k zz J9| in terms of the present model, the nodes in the B-phase 
have to be aligned parallel to [0 1]. 

The second objective of this work is to study and make predictions for the nodal exci- 
tations in the vortex state of the proposed model. Simple expressions will be derived for 
the quasiparticle density of states as well as for the angle dependent magnetothermal con- 
ductivity. The thermal conductivity obtained from this model describes the experimental 



data 



well, whereas the previously proposed order parameters do not. 



I. NODAL SUPERCONDUCTIVITY IN PrOs 4 Sbi 2 



Following Ref. [6[, we consider energy gap functions A(k) with point nodes at [0 1 0] 
and [0-10] in the B-phase, and with point nodes at [1 0], [0 1 0], [-1 0], and [0 -1 0] in 
the A-phase. Furthermore, in order to consistently describe the thermal conductivity data 
of Ref. 9( one also needs nodes at [0 1] and [00-1]. These constraints suggest 
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with e^ 1 = &2 ± iks etc. for the A-phase, and with 



(1) 



A B (k) = dAe 



±14 



1-kl 



(2) 



for the B-phase. We note that the p+h order parameter for the A phase satisfies the cu- 
bic symmetry. Also, for the B-phase the symmetry axis is rotated parallel to the crys- 
tal c-axis. Hence, the angular part, / = 1 — fc|, of the B-phase is the same as for 
the previously considered s+g superconductor, whereas the angular part of the A-phase, 
/ = | ( 1 — k\ — k\ — k\ J , is different. The magnitude of these superconducting order pa- 
rameters are shown in Fig. 1. Since both of these triplet order parameters break chiral 
symmetry, the ground state of the A-phase is sixfold degenerate and the ground state of the 
B-phase is twofold degenerate. 





FIG. 1: |A(k)| of the proposed p+h-wave superconducting order parameter in the A-phase (left) 
and in the B-phase (right) of PrOs4Sbi2- 



II. QUASIPARTICLE SPECTRUM 



In the absence of an external magnetic field, the quasiparticle density of states is given 



by 



g(E) = \x\Re{ 



v^H/T 



:)■ 



(3) 



where x = E/A, f 



1 — kf — k\ — k\ ) for the A-phase, and / = 1 — k\ for the B- 



phase. Here (...) = (47r) _ l f dQ... denotes the angular average. These quasiparticle densities 
of states are evaluated numerically and shown in Fig. 2. In particular for the low-energy 
limit |x| -C 1 we find g(E) ~ 7r|a;|/4 for the A-phase and g(E) w 7t*|x|/8 for the B-phase. 
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FIG. 2: Quasiparticle density of states in a p+h-wave superconductor, (a) A-phase, and (b) 
B-phase. 



In the vortex state the quasiparticle density of states at E = is given by ir(\v- g|)/(4A) 
(A-phase) and ir{\v-q\)/(8A) (B-phase), where v-q is the Doppler shift. J14j,[l5( Note that in 
this case (...) denotes the average over the unit cell of the vortex lattice and over the nodal 
lines and points on the Fermi surface. For the field dependence we find 



- lvVeH 

9a(H) = 2—^tAi 

-± lvVeH 

9b{H) = 1^— j b( 



(A — phase) 
(B — phase) 



(4) 
(5) 



with 

I A (0, 0) = sin 9 + (1 - cos 2 9 cos 2 0) 1/2 + (1 - cos 2 9 sin 2 0) 1/2 , (6) 

I B (M) = sin#. (7) 

For the B-phase it was assumed here that the point nodes are aligned parallel to [001]. 

The corresponding specific heat, the spin susceptibility, and the superfluid density in the 
clean limit and at ultra-low temperatures are given by 16[ 

C s /( lN T) = g{H), (8) 

Xs/xn = g(H), (9) 

p s (H)/p s (0) = 1 - g(H) (A - phase), (10) 

p 8 (H)/p s (p) = 1 - 3g(H) (B - phase). (11) 

The superfluid density in the A-phase is isotropic and given by Eq. 10. However, in the 
B-phase Eq. 11 is only valid if the supercurrent flows parallel to the z axis (parallel to the 
nodal direction). Therefore it will be of great interest to study the Knight shift in the vortex 
state. 

III. ANGLE DEPENDENT MAGNETOTHERMAL CONDUCTIVITY 

In order to analyze the thermal conductivity it is necessary to consider the effects of 
impurity scattering. Unlike for the s+g order parameter, the effects of disorder in the p+h- 
wave superconductor are more conventional. [l7|, [l8( Here we consider impurity scattering in 
the unitary limit. 

Let us first focus on the self-consistent equation for impurity scattering, given by iCoA = 
lim^^o^! where Cj is the disorder-renormalized frequency [151] For the A-phase, this leads to 

c ^{ c ^{-MTw) ) + {xt ^^) 1 • (12) 

where V is the quasiparticle scattering rate in the normal state, and x = \v ■ <f|/A. For the 
B-phase the prefactor 2 on the right-hand side of the above equation is replaced by 4. In 
the superclean limit (x) ^> Cq Eq. 12 leads to 

4r 

Coa = —r^y 1 (A-phase), (13) 

Cob = ^(x)- 1 (B-phase) (14) 

7TA 



On the other hand, in the clean limit we arrive at 

2 2F (r 2 ) 

Cl A In — = - - Y' ( A - phase ) ( 15 ) 

Cls In — = - - i-2, (B - phase). (16) 

In the absence of a magnetic field, the thermal conductivity exhibits universal heat 
conduction |l9l. |20j with k 00 /T = (7r 2 n)/(12mA) in the A-phase. In the B-phase there is 
universal heat conduction kqo/T = (ir 2 n)/(8mA) only when the point nodes are aligned 
parallel to the heat current J q . Here /too is the thermal conductivity in the limits T — > 
and T — > 0. Hence the experimental data of Ref. [9] are consistent with p+h-wave super- 
conductivity in PrOs 4 Sbi 2 , but inconsistent with a s+g-wave order parameter. 

In the vortex state, the thermal conductivity in the superclean limit is given by 

K zz v 2 (eH) 



k n 8A 

k zz _ 3v 2 (eH) 



., sin 2 9 (A-phase), (17) 

sin 2 9 (B-phase), (18) 



k n 64A 2 

where J q jj z. Here the thermal conductivity in the normal state is given by kn = 
(7i 2 nT) J (3mT) . For the B-phase we have assumed that the nodes are aligned to [0 1]. 
When the nodes are along [1 0] or [0 1 0], k zz is at least smaller by a factor of 10~50. 
Note that k zz ~ H sin 2 9 for both the A-phase and the B-phase. Therefore, the observed 
H-linear thermal conductivity k zz at T < 0.3f^ 9] follows from Eqs. 17 and 18. This implies 
that the crystals used in these measurements are in the superclean limit at sufficiently low 
temperatures. 

In the clean limit the expressions for the thermal conductivity become 



£ = i+ ^^hWwJ=) (A - phase) ' (i9) 



— = 1 + T^' n (^) 3in2gln ( 7W—J <B-i.l'-")- >-m 

k 00 12TA \\TJ \vVeHsm9/ 

where k qo is the thermal conductivity in the limit of universal heat conduction. Hence, 
in the clean limit the field-dependent part of k zz is given by k zz ~ HF(9), where 
F(9) = sin 2 9 ln(C/ sin 9) and C = A/vy/eH. From these equations we observe that the 
^-dependence of the leading terms is the same for the A-phase and B-phase. In Fig. 3, we 



plot the measured ^-dependence of k zz in PrOs 4 Sbi 2 |9j at various applied magnetic fields 
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FIG. 3: Angle dependence of the thermal conductivity in an applied magnetic field. The symbols 
are experimental data from Ref. |9|, and the solid lines are fits to HF(9), as discussed in the text. 

and azimuthal angles <fi. Fits of this data to HF{9) are shown as solid lines, suggesting that 
these samples are in the clean limit. In these fits C = 5 and C = 3 are obtained for H = 0.5T 
and H = 1.2T respectively. Thus, using the weak-coupling theory gaps A^ = 4.2K and 
As = 3.5K for the A- and B-phase, we deduce v = 0.96 x 10 7 cm/ sec and T ~ 0.1K. These 



values are reasonable. 



221 ] indicating that the quasiclassical approximation is reliable. 



IV. CONCLUDING REMARKS 



In this work, we have proposed a spin-triplet p+h-wave order parameter to account 
for observed features in the superconducting phases of PrOs4Sbi2- This model describes 
well the angle dependent thermal conductivity data by Izawa et al., Refs. [4] and J9(. 
In order to be fully consistent, we have discovered that the nodal directions of A(k) in 
the B-phase have to be aligned parallel to the external magnetic field in the field-cooled 
configuration. This triplet superconductivity in PrOs4Sbi2 is not surprising since many 
other heavy-fermion superconductors appear to have spin-triplet order parameters, including 
UPt 3 , UBei3, URu 2 Si 2 , and UNi 2 Al 3 . |2l( The interesting dependence of the nodal points 
in A(k) on the external magnetic field deserves further study 
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